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ABSTRACT 

We revisit the flare on 1992 January 13, which is now universally termed 
the "Masuda flare" . The revisit is motivated not only by its uniqueness despite 
accumulating observations of hard X-ray coronal emission, but also by the im- 
provement of Yohkoh hard X-ray imaging, which was achieved after the intensive 
investigations on this celebrated event. Through an uncertainty analysis, we 
show that the hard X-ray coronal source is located much closer to the soft X-ray 
loop in the re-calibrated HXT images than in the original ones. Specifically, the 
centroid of the Ml-band (23-33 keV) coronal source is above the brightest pixel 
of the SXT loop by ~ 5000 ± 1000 km (~ 9600 km in the original data); and 
above the apex of the 30% brightness contour of the SXT loop by ~ 2000 ± 1000 
km (~ 7000 km in the original data). We suggest that this change may naturally 
account for the fact that the spectrum of the coronal emission was reported to 
be extremely hard below ~20 keV in the pre-calibration investigations, whereas 
it has been considerably softer in the literature since Sato's re-calibration circa 
1999. Still, the coronal spectrum is flatter at lower energies than at higher en- 
ergies, owing to the lack of a similar source in the L-band (14-23 keV), which 
remains a puzzle. 

Subject headings: Sun: flares — Sun: X-rays, gamma rays — Sun: Corona 



1. Introduction 



As one of the most remarkable discoveries from the Yohkoh rn ission (jOgawara et al. 



199ll ). the flare occurring on 1992 January 13 (IMasuda et al.l 11994 the "Masuda flare") 
revealed for the first time a bright hard X-ray source in the corona, simultaneously with 
conjugate footpoint sources, which helped shape the modern vision of solar flares. Most 
importantly, it has corroborated the paradigm, long conjectured by theorists, that particles 
are accelerated high in the corona a nd stream down to the chromosphere along a coro- 



nal magnetic loop (see the review by lAschwanderu 120021 . and references therein). In spite 
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of accumulating observations of hard X-ray co ronal emission since then, the original Ma- 
suda flare has remained unique, as reviewed by Krucker et al. ( 2008 1*1. in that (a) the hard 
X-ray coronal source is located about 7000 km (~10") above the apex of the soft X-ray 
loop; and that (b) the hard X-ray coronal source is surprisingly weak in the 14-23 keV 
range. Its spectrum is t herefore extremely flattene d below ~20 keV, excluding a plausible 
thermal interpretation ( lAlexander fc Metcalj 119971 ). These features , however, are in con- 



trast to the coronal sources reported 



Battadia fc Benzl 12007 



Krucker fc Lin 



ater o n (e.g., 



Tomczak 2001: Petrosian et al. 2002 



20081 ). The survey on partially occulted flares ob- 



served by RHESSI ( Krucker fc Lin 20081 ). in which the footpoints are occulted by the limb 



but the thermal loop is visible, show that most hard X-ray coronal emission during the flare 
impulsive phase is only slightly above (< 6") the thermal loop, with many filled-loop events 
even displaying co-spatial non-thermal emission, and that the power-law index is between 4 
and 7, much softer than that of comparable on-disk flares. 

On the other hand, although the Masuda source has been described as being distinctly 
above the loop top, more recently, "above-the-loop-top" s ources have been pr i marily, if 



not exclusively, found in a double coronal source morphologv (ISui fc Holmanll2003l : ISui et al. 



20041 : iPick et al.ll2005l : IVeronig et al.ll2006l : iLi fc Gan]l2007l : iLiu et al.ll2008[ ) : the lower coronal 
source is located at the thermal loop top, while the upper sour ce is located as far as 30" 
above it, moving upward at a speed as high as 300 km s~^ (e.g., ISui fc Holmanll2003l ). The 



two sources "mirror" each other with respect to a presumed X-point reconnection site, in 
that higher energy emission comes fr om lower altitudes for the upper sou rce, while the lower 



source exhibits a reversed order (e.g.. lSui fc Holmarul2003l : iLiu et al.ll2008l ). This morphology 



is obviously different from that of the Masuda flare. 



Meanwhile, ISato et al.l (119991 ) accurately estimated the instrumental response functions 
through self-calibration using solar flares as calibration sources, and improved the Maximum 
Entropy Method (MEM) image-synthesis algorithm. Accordingly, the HXT images of the 
Masuda flare display somewhat subtle but significant changes, as compared with the pre- 
calibration publications (see Section 2). The spectral characteristics of the coronal source, 
as reported in the post-calibratio n investigations using the MEM algorithm, also change 
substantially ( iMasuda et al.l l2000l : iPetrosian et al.l l2002l also see Section 3). However, a 
mind-set has been developed, and these new developments are unfortunately left unaddressed 
or unnoticed in the solar community. 



'Their Figure 2 sho ws the Masuda flare with the new calibration (Sato et al. 1999I) . but their review 



is based on the work by iMasuda et all (jl994l ) and lAlexander fc MetcaHl (jl997( ) , therefore not reflecting the 
changes introduced by the re-calibration. See also §3.1. 
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Intending to resolve these i nconsistencies in th e hterature, we re-analyze the Masuda 
flare with both MEM and Pixon jMetcalf et al.lll996l : [Alexander fc Metcalj[l997n algorithms 
in this paper. In Section 2, we check the changes of the hard X-ray source locations, as 
compared with the original data. In Section 3, we review the historical investigations and 
interpretations (§3.1), and then redo the imaging spectroscopic analysis (§3.2). In Section 
4, we discuss the signiflcance of the geometrical changes (§4.1), interpret the spectral results 
from both the thermal (§4.2) and non-thermal (§4.3) viewpoints, and a simplifled picture is 
proposed to explain the morphological, as well as spectral, evolution of the coronal source 
(§4.3). 



2. Geometry 

Figure [1] shows the GOES soft X-ray fluxes (the top panel) and the HXT count rates 
(the bottom panel). The hard X-ray emission occurs in the impulsive phase and consists 
of a single peak that lasts for about 2 min. Three time intervals are chosen for image 
reconstruction, namely, 17:26:50-17:27:35 UT (slow rise phase), 17:27:35-17:28:15 UT0 (fast 
rise phase), and 17:18:15-17:29:00 (decay phase), which are marked by dotted lines. The 
Hard X-ray Telescope (HXT) provides four energy channels, L (14-23 keV), Ml (23-33 keV), 
M2 (33-53 keV) and H (53-93 keV). The coronal source is visible preferentially in the Ml- 
and M2-band. HXT images in the lower three channels are shown in Figure |2] as contours 
in the same color code as in the bottom panel of Figure [H Soft X-ray images obtained by 
the Soft X-ray Telescope (SXT) are displayed in grey colors. The three columns in Figure [2] 



Hhc same time interval used by Masuda et al. ( 2000[ ) 



Table 1. Flux Uncertainty 



Time Interval (UT) 


Algorithm 


L 


CS 


Ml 
NFP 


SFP 


cs 


M2 
NFP 


SFP 


17:26:50-17:27:35 


MEM 


1.7% 


8.1% 


5.6% 


6.7% 


13.2% 


14.2% 


11.5% 


Pixon 


1.7% 


5.7% 


7.8% 


7.4% 


14.5% 


15.6% 


9.7% 


17:27:35-17:28:15 


MEM 
Pixon 


1.3% 
1.5% 


7.0% 
7.7% 


3.8% 
4.3% 


3.5% 
4.6% 


12.4% 
14.6% 


5.0% 
3.3% 


6.0% 
6.7% 


17:28:15-17:29:00 


MEM 
Pixon 


1.2% 
1.0% 


6.0% 

8.6% 


2.8% 
3.6% 


3.8% 
5.0% 


13.2% 


2.9% 
3.6% 


3.1% 
5.4% 



Note. — CS stands for coronal source, NFP for the northern footpoint, and SFP for the 
southern footpoint. For the L-band, the flux is integrated over the whole loop. 
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Fig. 1. — Flare lightcurve in soft X-rays (top) and hard X-rays (bottom). Dotted lines mark 
three time intervals for the HXT image synthesis. The inset in the bottom panel shows one 
of the earliest SXT images for the Masuda flare, which was taken at the onset of the flare 
impulsive phase. 
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Fig. 2. — Masuda flare observed with the re-calibrated HXT images. Hard X-ray contours 
are shown in the same color code as in Figure [H with the SXT images as the background. 
Contour levels are 10.0, 20.0, 40.0, and 80.0% of the maximum brightness in each individual 
energy band. The three columns corresponds to the three time intervals indicated in the 
bottom panel Figure [H Images in the top (bottom) row are reconstructed with the MEM 
(Pixon) algorithm. 
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Table 2. Geometry of the Masuda Flare 



Interval 


Algorithm 


HsxT^ 
(10^ km) 


H'sxt'^ 
(10^ km) 


Hl" 
(10^ km) 


(10^ km) 


(10^ km) 


Dppf 
(10^ km) 


(10^ km) 


(10^ km) 


Slow Rise 


MEM (sim) 
Pixon (sim) 
Pixon (err) 


136±9 
132±9 
129±9 


171±9 
167±9 
164±9 


138±4 
131±5 
130±3 


183±6 
177±7 
182±6 


212±9 
209±8 
203±5 


173±4 
179±5 
182±7 


47±11 
45±12 
53±11 


12±11 
10±12 
18±11 


Fast Rise 


MEM 
MEM (sim) 
Pixon (sim) 
Pixon (err) 


125 

143±9 
143±9 
142±9 


178±9 
178±9 
177±9 


148 

143±2 
145±6 
136±6 


221 

194±4 
193±5 
192±3 


213±9 
211±9 
216±11 


157 

173±3 
169±3 
172±2 


96 

51±10 
50±10 
50±9 


70 
16±10 
15±10 
15±9 


Decay 


MEM (sim) 
Pixon (sim) 
Pixon (err) 


140±9 
141±9 
141±9 


175±9 
175±9 
176±9 


149±2 
139±4 
138±8 


201±8 
197±8 
198±6 


222±7 


173±2 
174±2 
175±3 


60±12 
56±12 
57±11 


26±12 
22±12 
22±11 



Note. — In the brackets of the 2nd column, "sim" and "err" indicates that the uncertainties are derived from the simulated 
HXT images, and from the Pixon error map, respectively. Fo r comparison, the measurements for the origi nal Masuda flare from 
17:28:04 to 17:28:40 UT made by I Aschwanden et al] |l9963) are listed in bold face; the measurement bv lMasuda et al.l l ll994l ) 
is given in the last column in slanted face. 

''The height of the brightness maximum of the SXT loop, which is represented by the its distance to the midpoint of the 
centroid positions of the Ml-band conjugate footpoints. A similar approach is taken for other columns. Also see Figure [3] 

''The height of the apex of the SXT loop. 

'^The height of the centroid position of the enhanced L-band loop top. 
''The height of the centroid position of the Ml-band coronal source. 
''The height of the centroid position of the M2-band coronal source. 
'The distance between the conjugate Ml-band footpoints. 

^The height difference between the Ml-band coronal source and the brightness maximum of the SXT loop. 
''The height difference between the Ml-band coronal source and the apex of the SXT loop. 
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correspond to the three time intervals indicated in Figure [TJ Images in the top (bottom) row 
are reconstructed with the MEM (Pixon) algorithm. 

One can immedi ately see in Figu re |2] that the L-band loop is different in shape from 
the original one (cf., ISato et al.l Il999l ). and that the hard X-ray coronal source (Ml- and 
M2-bands), located only slightly above the SXT loop, is generally enclosed by the L-band 
loop at 10-20% brig htness contour level jjl. In contrast, in the original MEM image (e.g., left 
panel of Figure 3 in lAschwanden et al.l |l996aj) , the hard X-ray coronal source was spatially 
distinct not only from the SXT loop, but from the L-band loop as well. 



Further comparisons are made quantitatively. In a similar approach taken by lAschwanden et al. 



(Il996aj), we measure the centroid position of the loop-top source, and compute its height, 
in each HXT energy band (except the H-band), with r espect to the conjugate footpoints in 
the Ml-band. For SXT, following lAschwanden et al.l fll996a| ). the brightness maximum at 
the enhanced loop top is pinpointed, for which the measurement error is estimated to be 
hal f of the SXT CCD p ixel size (2.46", or ~1750 km). To compare with the result given 
by iMasuda et al.l (11994! ) that the Ml-band coronal source was above the apex of the SXT 
loop by ~7000 km, we draw a line connecting the midpoint of the Ml-band footpoints and 
the brightness maximum of the SXT loop, and choose the pixel where the line intersects the 
SXT 30% brightness contour (see Figure |3]) as a conservative proxy for the loop apex. 

To estimate the uncertainties of the hard X-ray source locations, we simulate the HXT 
images by adding random noise to the data counts collected by each of the 64 subcoUimators, 

Cobs; i-G-; 

IS + \/ Ci^ obs ^ "^ii 



obs 



1, 2, 



64, 



where is a pseudo-random number generated by the IDL procedure RANDOMN. A new image 
can be reconstructed based on Csim and the source centroids are measured. In practice, a box 
is specified to enclose each individual source in an image. A region is allowed to grow from 
the local maximum inside each box to include all connected neighbors whose values are above 
a given minimum (see the REGION_GROW function in IDL), which is chosen interactively so 
that the region takes the shape of an ellipse. The centroid position is determined by fitting a 
2D elliptical Gaussian equation to each individual region. By repeating the above procedure 
for many times (30 times in our work), we get a distribution of the source centroids, and 
their uncertainties can be approximated by standard deviation. The flux uncertainty of each 
individual source as shown in Figure |2] can also be obtained from the simulated images (see 



■fThe nominal dynamic range of HXT is 1:10. Due to the nonlinear nature of both MEM and Pixon 
algorithms, however, the image quality is not only dependent on the observation error/noise, but on the 
source distribution/complexity as well. 
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Table [T]) by integrating the flux inside a box region containing the source. 



For the Pixon algorithm, alternatively, an error map (l-cx error at each pixel) can be 
obtained for the reconstructed image. Thus, we can estimate the uncertainties of the source 
centroids by adding (substracting) the error map to (from) the reconstructed Pixon map, 
getting the "new" source locations, and then comparing with those obtained from the Pixon 
map alone. The larger of the differences is taken as the error bar. 

In Figure [31 we show the source centroids during the fast rise phase (17:27:35-17:28:15 
UT), with corresponding error bars estimated from simulated MEM (left panel) and Pixon 
(middle panel) images, as well as from Pixon error maps (right panel). Sources in different 
HXT energy bands are represented with the same color code as in Figure [H The brightness 
maximum of the SXT loop top is shown in orange, and its apex is denoted by an asterisk. 
One can see that the accuracy of the HXT loop-top centroid position is comparable to that 
of the brightest SXT pixel, despite that the FWHM of the finest HXT collimator is as large 
as 6000 km. The footpoint positions are even more accurate, with sub-arcsec precision being 
achieved. 

Various geometrical parameters can then be derived and their error bars are computed 
following the rules of the propagation of unce rtaint^yjl. The results are lis ted i n Table [21 For 



compa rison, corresponding numbers given by lAschwanden et al.l (Il996al ) and iMasuda et al. 



( 119941 ) are also listed in bold and slanted face, respectively. One obvious difference is the 
smaller separation of Ml-band coronal source from the brightness maximum of the SXT 
loop. The height difference, AH in Table [2 (5 000 ± 1000 km in Pixon a nd 5100 ± 1000 
km in MEM), is almost half of that measured by lAschwanden et al.l (Il996ai 9600 km). The 
height difference between the Ml-band coronal source and the apex of the SXT loop, AH' in 
Table [21 (1500 ± 1000 km in Pixon and 1600 ± 1000 km in MEM), is also much smaller than 
that given by iMasuda et al.l (Il994l. 7000kni ) . The new measurement is in agreement with 
RHESSI observations ( iKrucker Sz Linll2008l ) that typical impulsive phase coronal emission 
is only slightly above the thermal loop (< 6"), which, however, is measured by the radial 
dif ference of the cente r of ma ss locations. This is slightly different from the approach taken 
by lAschwanden et al.l (Il996al ) and this paper. 



Accordingly, the Ml-band coronal source is also less separated from the L-band loop top 
in terms of centroid positions. Thus, any reasonably chosen integration region containing the 
Ml-band source in imaging spectroscopic analysis will inevitably include much more L-band 



''A note of caveat sho uld be made that th e SXT-HXT co-aUgnment error is not included. However, it is 
generally better than 1" (jMasuda et aLlll995[ ). less than the error of locating the brightest SXT pixel (^900 
km). 
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photon counts than it would with the original HXT images, therefore leading to a much 
larger L- and Ml-band ratio, as we shall demonstrate in Section 3. 



3. Spectroscopy 



3.1. Historical Investigations and Interpretations 

Limited spectral information from HXT images can be derived by calculating flux ra- 
tios between adjacent energy bands. We denote the photon spectral index derived from L- 
and Ml-band as 71, and that from Ml- and M2-band as 72, with the superscripts, CS and 
FP, standing for the coronal source and the footpoints, respectively. Corresponding electron 
indices are denoted by 5 in a usual fashion. From Table P one can see that the coronal spec- 



trum is considerably softer in the pos t- calibration era (iMasuda et al 



2OOOI: IPetrosian et ah 



2002) than in the pre-calib ration era ( Masuda et al. 19941 19951 : Alexander fc Metcall 



1997 



Metcalf fc Alexanderl I1999I ). Different decisions made by data analyzers in the same era, 
e.g., different integration region and time interval for data accumulation, appear to result in 
only minimal differences, however. 



In the pre-calibration era, the verdict for the Ma suda flare is that the coronal emis- 
sion was nontherma l (e.g.. [Alexander fc Metcalj Il997l ). The thermal interpretation (e.g., 
Masuda et al.l 119941 ) was largely rejected due to the following three major arguments: (a) 



the thermal bremsstrahlung requires an extremely hot plasma (up to T > 200 MK), for 
which there was no evidence from other instruments, e.g. the impulsive pha se BCS Fe XXV 
spectra, assumed from the loop top, give a source temperature of ~ 30 MK (lFletcherlll999l ): 
(b) the temperature derived from the HXT lo wer channels is much higher (> 200 MK) than 
that from the higher channels (125±15 MK; [Alexander fc Metcalj 119971 ): and (c) the time 
variability of the coronal source wa s too rapid (similar to that of the footpoints) to be con- 
sistent with thermal cooling times ( iHudson fc RyanI 119951 ). We will check these arguments 
against the re-calibrated data in §4.2. 

On the other hand, it is very difficult to understand the extremely hard spectrum of the 
coronal source, which had a power-law index of ~2 from the L- and Ml-band ratio. Though 
it may indicate a possible turnover in a nonthermal electron spectrum, such a fiat spectral 
index is close to the hardest theoretically possible brem sstrahlung spectrum. Never theless, 
putting in the perspective of the thick-thin target model (jWheatland &: Melroselll995l ). which 



argues that a coronal source acts as a thick target {S = 7f -|- 1) for electrons at lower 
energies and thin target (5*"^ = 7^^ — 1) for higher energies, one can derive that the spectral 
index of the injecting electrons for the coronal source, S^^ ~ 3, while that for the footpoint 
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Fig. 3. — Positions of HXT and SXT sources during the fast rise phase (17:27:35-17:28:15 
UT). Error bars of HXT sources are estimated (see the text) from simulated MEM (left panel) 
and Pixon (middle panel) images, as well as from the Pixon error map (right panel). Source 
centroids in different HXT energy bands are shown with the same color code as in Figure [1] 
The brightness maximum of the SXT loop is shown in orange, whose error of measurement 
is half of the SXT CCD pixel size (1.23")- The SXT loop observed at 17:27:43 UT is shown 
in grey colors, and overlaid by the contour at the 30% of its brightness maximum. The apex 
of the SXT loop is denoted by an asterisk (see the text). 



Table 3. Historical Numbers 



Time Interval (UT) 


Algorithm 


L/Ml 


M1/M2 


7 

L/Ml 


FP 

M1/M2 


tcs (mK) 
L/Ml M1/M2 


17:26:52-17:27:39^ 


MEM 


2.6 


4.1 






200 


130 


17:27:40-17:28:20^ 


Pixon 


2.2±0.f 


3 4.1±0.2 


3.3±0.1 


3.5±0.1 


215±60 


120±15 


17:27:35-17:28:15^ 


MEM 


4.0 


5.5 


3.2 


3.8 


100 


90 


4 


MEM 


4.2 


6.8 


3.1/2.5 


3.9/3.7 







1 From iMasuda et al] lll994 Il995l) 

^From [Alexander fc Metcalj lll997l l and lMetcalf fc Alexander! |l99i) 
^From lMasuda et al.l 1 I2OOOI) 

*From IPetrosian et al. I l l2002h . who specified only a time instant, 17:28:04 UT, but not the time 
interval for data accumulation. Note also that is derived for each individual footpoint. 
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emission, S^^ ~ 4.5. Note, however, that the thick-thin target model assumes that the 
same electron population is responsible for both the coronal and the footpoint emission. A 
promising explanation is electron tra pping in the corona, with energy losses dominated by 
collisions ( lAlexander 8z Metcalj 119971 ). Since the collisional deflection time of an electron of 
energy E, 

E \^^^ / n \ -1 / 20 
t{E) ^ ^ " ^ ^ 



0.95 



1 keV 



108 



cm 



In A 



where Ue is the electron density and In A the Coulomb logarithm (~ 20 under solar con- 
ditions), the longer life time of more energetic electrons would result in the erosion of the 
low-energy regime of the velocity distribution function, therefore the progressive spectral 
hardening, on timescales in the order of 10-100 s for a coronal density of order 10^ cm~^. 

The coronal trapping of energetic electrons predicts that the injecting spectrum for the 
coronal source is hardening for 1.5 powers relative to the corresponding footpoint emission 
in the weak diffusion limit (see §4.3). It also leads to the progressive delay of hard X- 
ray p eaks with increasing ener gy, due to collisional precipitation of the trapped electrons 



(e.g.. lAschwanden et al.lll996bl ). Based on time-of- flight measurements, lAschwanden et al. 



(jl996bl ) found 5 hard X-ray coronal sources that can be expla ined by a coronal trap in th e 
weak diffusion limit. The same 5 events were also studied by iMetcalf fc Alexanderl ( 19991 ). 



who were looking for the spectral evidence of trapping. Two of them, including the Masuda 
flare, were found to be consistent with impulsive phase coronal trapping in the thick-thin 
target scenario. 

In the post- calibration era , in contrast, the spectrurn of the Masuda loop-top source 



looks rather soft flMasuda et al.ll2000l : iPetrosian et al.l 12002), with its spectral in dices falling 



in the typical ra nge of recent RHESS I observations (4-7; iKrucker fc Linll2008l ). Using the 
MEM algorithm, iMasuda et al.l (120001 ) reported that the loop-top spectrum can be well fit 
by emission from an isothermal plasma of about 100 MK. In §3.2, we present a re-analysis 
of the Masuda flare with both MEM and Pixon algorithms, and in Section 4, discuss the 
implications of the changes introduced by the re-calibration. 



3.2. Re- Analysis 

The region growing method (see Section 2) is also used to specify the integration region 
for each individual hard X-ray source in imaging spectroscopic analysis. This is a convenient 
method to "extract" weak sources, and the same integration region can be recovered as long 
as the same minimum threshold is specified. Integration regions chosen for the fast rise phase 
of the flare are specified as an example in Figures H] and O In Figure H] the integration regions 
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Table 4. Spectra Characteristics of the Masuda flare with the Integration Regions 
Covering both the Ml- and M2-band Emission 



Time Interval (UT) 


Algorithm 


7 

L/Ml 


cs 

M1/M2 


7 

L/Ml 


FP 

M1/M2 


L/Ml 


tcs (mK) 
M1/M2 


L/M1/M2 


7^ 


17:26:50-17:27:35 


MEM 
Pixon 


4.5±0.7 
4.6±0.3 


5.3±1.3 
5.0±0.4 


3.3±0.4 
3.2±0.3 


4.6±0.5 
4.0±0.2 


80± 18 
78 ±7 


92 ± 34 
101 ± 14 


84± 15 
86 ±6 


4.03±0.06 


17:27:35-17:28:15^ 


MEM 
Pixon 


4.5±0.4 
4.7±0.2 


5.9±1.1 
5.8±0.5 


2.5±0.2 
1.9±0.2 


3.7±0.2 
3.3±0.1 


79± 12 
73 ±7 


78 ± 21 
80± 10 


79 ±9 
76 ± 5 


3.69±0.02 


17:28:15-17:29:00 


MEM 
Pixon 


5.8±0.7 
6.6±0.4 


5.6±1.8 
10.0±5.1 


1.8±0.2 
1.3±0.5 


3.5±0.1 
3.1±0.1 


54 ±9 
45 ±4 


84 ±41 
37 ±26 


57 ±9 
45 ±4 


3.36±0.02 



^Spectral indices derived from single power law fits of the corresponding Yohkoh Wide Band Spectrometer (WBS) spectra. 
■^Integration regions are specified in Figure |4l 



Table 5. Spectra Characteristics of the Masuda flare with the Integration Regions 

Covering the M2-band Emission 



Time Interval (UT) 


Algorithm 


7 

L/Ml 


CS 

M1/M2 


7 

L/Ml 


FP 

M1/M2 


L/Ml 


T^s (MK) 
M1/M2 


L/M1/M2 


17:26:50-17:27:35 


MEM 
Pixon 


4.3±0.8 
4.4±0.3 


4.8±1.4 
4.8±0.4 


3.3±0.4 
3.2±0.3 


4.6±0.5 
4.0±0.2 


84 ±24 
81 ±8 


105 ± 46 
104 ± 15 


91 ± 19 
89 ±6 


17:27:35-17:28:15^ 


MEM 
Pixon 


3.0±0.7 
3.4±0.5 


5.4±1.3 
5.1±0.5 


2.3±0.2 
1.8±0.2 


3.6±0.2 
3.3±0.1 


152 ± 61 
124 ± 30 


88 ±31 
97± 13 


117 ± 23 
106 ±9 



'^Integration regions are specified in Figure [5] 
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Fig. 4. — Integration regions for the fast rise phase. 1-sigma error maps of the corresponding 
HXT images are shown in 2nd and 4th rows. For MEM (Pixon), the integration regions cover 
the union of the 10% contour of the maximum Ml- and M2-band emission for the footpoints, 
and the union of the 10% contour of the maximum Ml-band emission and the 6% (10%) 
contour of the maximum M2-band emission for the loop-top source. 
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cover emission from both the Ml- and M2-band, and the corresponding results are hsted in 
Table HI in Figure |5l integration is carried out over the regions with significant M2-band 
emission only, and the corresponding results are listed in Table |3 Table O does not include 
the decay phase, during which there is no visible M2-band coro nal emission. From both 
tables, one can s ee tha t our results are basically in agreement with iMasuda et al.l ( l2000l ) and 



Petrosian et al.l (120021 ) . In addition, one can see that the photon spectra of the footpoint 



emission are hardening throughout the hard X-ray burst. In contrast, the spectra of the 
coronal emission are successively softening, which we shall discuss in §4.3. 



4. Discussion and Conclusion 



4.1. Impact of the Geometrical Change 



There are three major difficulties in understanding the original Masuda fiare, viz., 1) the 
lack of L-band emission at the coron al source and consequently the extremely flat sp ectrum 
of the coronal source below ~20 keV (IMasuda et al.ll 1994 [Alexander fc Metcalilll997l ): 2) the 
apparently low density at the location of the coronal source due to the lack of SXT emission; 
and 3) the lack of soft X-ray signature of hydrodynamic response at th e presumed fleld lines 
that connected the hard X-ray coronal source to the chromosphere (j Alexander fc Met call 
19971 ). Also, the soft X-ray loop barely grows with time, defying the theory of choromospheric 
evaporation. All three difficulties can be arguably attributed to the elevated hard X-ray 
coronal source in the original data. Thus, the reduction of the source height relative to the 
soft X-ray loop in the re-calibrated data has more or less relieved these difficulties. 

We will discuss the 3rd difficulty in §4.3. As for the first one, although the spectrum 
has become softer at lower energies with the re-calibrated data, the fact that the Masuda 
source is poorly imaged in the L-band remains a puzzle. It might be due to a turnover in 
non-thermal electron spectrum, or to instrumental effects, or simply to the poor resolution 
of HXT. In any case, the L- and Ml-band ratio reflects poorly the true spectral slope. 

The 2nd point deserves some further comments. The column depth required to stop 
electrons at 20 keV, the approximate break energy of the coronal spectrum, is N ~ 10^*^ 



cm ^, which can be estimated from the formula (e.g., IWheatland fc Melrosd 119951 ) 



E = {2KN) 



1/2 



where E is in erg, N is in cm~^, and K = 27re^lnA, with e the electron charge (4.8 x 10^^*^ 
esu). Given a source length of / 2^ 10^ cm, the density, rip, is as large as 10^^ cm""^. This, 
however, agrees approximately with lAschwanden et al.l (119971 ) who found the electron density 
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in the flare trap of tlie Masuda flare, — (0.21 ±0.04) x 10^^ cm ^, by fitting the hard X-ray 
energy-dependent time delays with the electron coUisional deflection time. The trap density is 



lower than the peak density in the SXT loop top, (0.96±0.20) x 10 cm (lAschwanden et al. 



19971 ) , but at least an order of magnitude larger than the de nsity at the original co ronal source 



location, e.g., 3 x 10^ cm ^, inferred from HXT data by iMasuda et al.l (119941 ). or 3 x 10 



cm ^, inferred from both HXT and SXT data by iTsuneta et al.l (119971 ). Hence the coronal 



trap must be located above the brightness maximum of the SXT loop, but lower in altitude 
than the original Masuda source. If the trap location is consistent with the current source 
location, the condition for the confinement of nonthermal electrons at the loop-top region is 
much less stringent than that imposed by the original data. For example, a 20 keV electron 
needs to bounce only 5 times for the thick target hypothesis to be effective, in contrast to 
more than 30 times of bo unce as required by the original data, given the column depth. 



N = nj^3x 10^^ cm-2 (iTsuneta et al.lll997n 



In view of the argument above, the thick-thin target model proposed by I Wheatland fc Melrose 



(11995! ) looks quite promising. It is noteworthy that the SXT loop is clearly visible at the 
flare onset (the bottom panel of Figured]), suggesting that the coronal atmosphere has al- 
ready been heated to a dense and hot state prior to impulsive electron acceleration. If the 
corona is dense enough, it can stop the less energetic electrons to produce a bright, thick- 
target coronal source in the L-band, and a much dimmer, thin-target one in higher-energy 
bands. RHESSI has observed coronal sources brighter than footpoints at energies up to 
50 keV (IVeronig fc Brownll2004l ). which take the shape of a loop and are explained by the 
thick-target model. We will test the re-calibrated data against the thick-thin target model 
(§4.3), but before that we first re-examine the thermal interpretation (§4.2). 



4.2. Thermal Interpretation 



From Table HI one can see that a consistent temperature can be derived from two 
different HXT chan n el ratios (L/Ml and M1/M 2), in contrast to previous publications 
(IMasuda et al.l Il994l : I Alexander fc Met calfll 19971 1. In Table [Sj temperatures derived from 
the lower HXT channel ratio are indeed higher than those from the higher channel ratio 
during the fast rise phase, probably suggesting the existence of a nonthermal component, 
despite that the two temperatures still fall within the same error range. 



We now check another major argument against the thermal interpretation, viz., the 
time variability of the coronal source is much shorter than the electron thermalization time 
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for Coulomb collisions (IHudson fc Rvanlll995h 



= 6 X 10" 



MK ) 



3/2 



10^ cm~3 



With np ^ 3 X 10^ cm ^ (assuming a unity filling factor) and T ~ 200 MK from lMasuda et al. 



(jl994j ). Tee — 60 s. This number has been commonly cited in the literat ure, out of the 



contex t of specific plasma parameters. With rig ^ 3 x 10^ cm~^ estimated by iTsuneta et al. 
( 119971 ) and T ^ 80 MK from our analysis, however, Tp,, is about 1.4 s. If the trap density 
(?7,e — 0.2 X 10^^ cm~^) inferred by lAschwanden et al.l fll997l ) is adopted, r^e could be as small 



as 0.2 s. Thus, the thermal interpretation cannot be ruled out unambiguously, based on the 
spectral information available. However, th e temperature is still much higher than that 
derived from BCS Fe XXV spectra (~ 30 MK: lFletcherlll999[ ). Moreover, a si ngle thermal fit 
can b e excluded for most hard X-ray coronal sources observed by RHESSI f lKrucker fc Lin 
20081 ). These observations cast doubt on the thermal interpretation of the Masuda fiare. 



4.3. Thick-Thin Target and Trapping 



Following the thick-thin target model (jWheatland fc Melroselll995l ). we assume that the 



coronal source acts as a thick target for less energetic electrons = + 1)? and a thin 
target for more energetic electrons (5^^ = 7^^ — 1). With low-energy electrons stopped 
in the loop-top region, the footpoint emission is not supposed to follow power laws at low 
energies, but a thick target is assumed at high energies as usual {5^^ = 7^^-1-1). The derived 
electron indices are given in Table El in which we only use the photon indices obtained by 
Pixon. One should trust (5^^ more than Sf^, since the latter is inevitably contaminated by 
thermal emission, especially for the integration region that covers both the Ml- and M2-band 
emission. The thick-thin target model agrees well with the observation during the fast rise 
phase of the hard X-ray burst, when 6^^, 62^, and 6^^ are approximately equal to each other 
(shown in bold face in Table E]). During the slow rise phase, however, is significantly 
smaller than S^^ . 

One difficulty in understanding the Masuda fiare (see §4.1) lies in the lack of signature 
of choromospheric evaporation. This can be understood if the less energetic electrons, which 
constitute the bulk energy of injecting electrons, are mostly stopped in the dense corona, 
as indicated by the existence of the soft X-ray loop at the fiare onset (Figured]). In that 
case, chromospheric evaporation would be strongly suppressed, due not only to reduced 
energy deposited in the chromosphere, but also to enhanced inertia of the overlying material 



( lEmslie et al.lll992[ ). Moreover, the fact that the electrons interact with hot coronal plasma 
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1000 



X (arcsecs) X (arcsecs) X (arcsecs) 

Fig. 5. — Integration regions for the fast rise phase. For MEM (Pixon), they cover the 6% 
(10%) contour of the maximum M2-band emission. 




Fig. 6. — Cartoon illustrates that the magnetic field strength at the hard X-ray loop-top 
region increase with time due to the field line shrinkage. The hard X-ray footpoint sources 
are indicated by blue elliptical shapes, and the hard X-ray coronal source is indicated by a 
red elliptical shape. Field lines are shown in solid line with arrows. The thin, orange loop 
indicates the 10% contour of the L-band emission, and the thick, yellow loop indicates the 
soft X-ray loop. 
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rather than a "cold target" reduces the energy-loss rate, hence fewer e lectrons are required 
at low energies, which helps to alleviate the electron number problem ( lEmslid l2003l ). 



The spectral characteristics revealed in Table E] seem to put the coronal trap in the weak 
pitch angle diffusion regime, in which the loss cone is empty, while the opposite case, that 
the loss cone is filled with scattered particles, is defined as the strong diffusion. Specifically, 
in a trap-and-precipitation scenario, the precipitation rate, z/, can be evaluated as. 



V cx 



n. 



weak diffusion limit, 
strong diffusion limit. 



where 6*0 is the loss-cone angle. If the data accumulation time is significantly longer than the 
trapping time (1/z^), which is often the case, the observed injecting electron spectrum for 
the coronal source in relation to that for the c orresponding footpoint emis sion would satisfy 
the following inequality (c.f., the Appendix in lMetcalf fc Alexanderlll999( ). 



-3/2<5^«-5^^<l/2, 
with the lower (higher) bound corresponding to the weak (strong) diffusion limit. 

The decrease in the additional hardening of ^""^ relative to with time, as indi- 
cated in Table O may suggest an enhanced scattering rate. Alternatively, this can be at- 
tribu ted to an interesting effect in the standard flare model, namely, field line shrinkage 



(e.g.. iForbes fc Actorull996l ). The cartoon in Figure E] shows that the magnetic field strength 
at the soft X-ray loop top would increase with time, as the newly reconnected, cusp-shaped 
field lines "shrink" and relax into a more potential state, piling up above the soft X-ray loop. 
As we know, the loss-cone angle, Qq^ is defined through the formula, 

sin2^o = 5/5o, 

where Bq and B are the magnetic field strength at the footpoint and at the coronal trap, 
respectively. The increase of the magnetic field strength at the coronal trap would result in 



Table 6. Power-Law Indices of Injecting Electrons Assuming the Thick-Thin Target Model 



Time Interval (UT) 


Integration Region 


<5p 


5P 


^FP 


17:26:50-17:27:35 


Ml & M2 


5.6±0.3 


4.0±0.4 


5.0±0.2 










M2 


5.4±0.3 


3.8±0.4 


5.0±0.2 




Ml & M2 


5.7±0.2 


4.8±0.5 


4.3±0.1 


17:27:35-17:28:15 


M2 








4.4±0.5 


4.1±0.5 


4.3±0.1 
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the enlargement of the loss cone, and therefore the enhancement of the precipitation rate, 
if we assume that electrons are continuously injecting into the trap. Electrons that are 
trapped before the field line shrinkage, however, would not be affected by the strengthening 
B, as their perpendicular momenta also increase due to the conservation of the transversal 
adiabatic invariant, p\/B. 

With the depletion of energetic electrons, the hard X-ray coronal emission would become 
weaker relative to the footpoints. Its centroid location would move to a slightly higher alti- 
tude relative to the soft X-ray loop, because electrons trapped at lower magnetic loops (with 
larger loss-cone angles) would be depleted from the coronal trap first. With the enhanced 
precipitation rate, low-energy electrons would still be stopped at the thick-thin target, but 
more high-energy electrons would precipitate towards the chromosphere. Thus, the footpoint 
spectrum would be successively hardening, while the coronal spectrum successive softening, 
as observed in the Masuda flare. 

To summarize, we demonstrate through an uncertainty analysis that the hard X-ray 
coronal source of the Masuda flare is located much closer to the soft X-ray loop in the re- 
calibrated HXT images than in the original ones. This may account for the change of the 
spectral characteristics of the coronal source derived from imaging spectroscopic analysis. 

We are grateful to the referee, Sam Krucker, for crucial help on the uncertainty analysis, 
as well as many thoughtful comments and suggestions, which substantially improve the 
quality this paper. RL thanks Wei Liu for helpful comments. This work was supported by 
NASA grant NNX08-AJ23G and NNX08-AQ90G, and by NSF grant ATM-0849453. 
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